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ABSTRACT 


Geostron’iic  relative  vorticit;/  is  confuted  along  three  latitude 
circles  generally  encompassing  the  United  States,  using  a  Fourier  analy¬ 
sis  of  the  500-rab  height  field  for  the  period  3-6  November  19^1 . 

A  statistical  comparison  of  vorticity  obtained  in  this  manner  with 
it*s  square  grid  computation  is  made. 

The  writers  are  deeply  indebted  to  Professor  F.  L,  Martin  of  the 
United  States  Naval  Postgraduate  School  for  his  suggestion  of  the  topic 
and  his  continued  help  throughout  the  investigation  and  during  the 
preparation  of  this  paper* 
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Symbol  Definition 


% 

vertical  component  of  geostrophic  relative  vorticity 

g 

acceleration  due  to  gravity 

f 

coriolis  parameter 

z 

contou.r  height 

a 

radius  of  earth 

0 

angular  velocity  of  earth 

latitude 

a  r 

difference  in  latitude 

* 

longitude 

4  JGL  - t;  f/  i 

'IT 

zonally  averaged  contour  height 

Z  " 

perturbation  of  zonal  contour  height 

z 

space -mean  contour  height 

d 

H 

grid  mesh  =  381  km 

wave  number 

latitude  circles  numbered  1  through  27  grid  lengths  from  the 
north  pole 

Cn 

amplitude  of  nth  wave 

©n 

phase  angle  of  nth  wave 

m 

map  factor 

r 

radius  of  earth  at  latitude  <b  ,  ^  -  d  Cj;  ‘ 

J 

F/V//F 

vertical  component  of  geostrophic  relative  vorticity  as  com¬ 
puted  by  Fleet  Numerical  Weather  Facility 

F 

F  test  for  equality  of  variances 

p 

number  of  predictors  (independent  variables ) 

ug 

trornic;-liy  -computed  zonal  component  cf  the  wind 

r 

i— o.l  v r. .t  .ed  geostrophically 

iii  $ 


Symbol 


Definition 


RKS 

R 


root  mean  square 

multiple  correlation  coefficient 
number  of  cases  in  sample 
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1 .  Introduction. 


The  advent  of  the  electronic  computer  coupled  with  increased  data 
coverage  has  removed  the  primary  obstacles  to  the  use  of  harmonic  analy¬ 
sis  as  a  useful  instrument  in  atmospheric  research . 

Use  is  made  of  a  Fourier  analysis  of  the  500-mb  height  field.  A 
computer  program  devised  by  G.  Arnason  and  M.  Reese,  and  made  available 
by  Fleet  Numerical  Weather  Facility"'  was  used  to  obtain  the  harmonic 
analysis  around  latitude  circles  at  distances  of  1  through  27  grid  lengths 
d  =  381  km  from  the  pole.  An  equation  for  relative  vorticity  using  the 
results  of  the  spherical  harmonic  analysis  is  derived  and  tested  statisti¬ 
cally  against  the  more  usual  numerical  computation  of  geoctropbic  relative 
vorticity  made  using  a  square  grid. 

The  problem  undertaken  was  to  compute  relative  vorticity  at  500  mb 
using  the  peostrophic  vorticity  equation  in  its  spherical  coordinate  form 
[1,  P-  358], 
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In  equation  (l  ),  z  was  replaced  by  its  expression  in  terms  of  a  Fourier 
analysis 

7  ~  i  t  £[_C,,  -e\)'J  -  f  A-  /*  (2) 

n;( 

Here  z  is  the  zonally  averaged  contour  height,  whereas  z*‘  is  the  perturba¬ 
tion  quantity  which  remains  when  the  z  field  is  subtracted  from  the  actual 
height  field.  The  computer  printout  included  maps  of  each  of  these  three 
fields  at  500  mb  for  the  period  under  investigation,  3  through  6  November 
1961.  Geostrophic  relative  vorticity  was  computed  by  the  above  mentioned 
technique  for  17  longitudes  along  each  of  three  latitude  circles,  generally 
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Henceforth  referred  to  as  FNVJF 
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encompassing  the  north-south  extent  of  the  United  States.  This  gave  a 
total  snmnle  size  of  204  computed  relative  vorticities  which  were  then 
correlated  with  square-grid  computations  of  relative  vorticities  obtained 
from  FNUF  charts  of  absolute  vorticity  for  exactly  the  same  set  cf  points 
and  times. 

2.  The  geostrophic  relative  vorticity  in  terms  of  a  harmonically  analyzed 
contour  field. 

The  adequacy  of  the  Fourier  representation  of  meteorological  para¬ 
meters  will  not  be  discussed  here;  the  interested  reader  is  referred  to 

L2J. 


It  was  previously  mentioned  that  use  was  made  of  the  FNWF  computer 
program  for  a  Fourier  analysis  of  the  500  mb  contour  fiold.  A  brief  de¬ 
scription  of  the  theory  and  numerical  evaluation  of  the  Fourier  coeffi¬ 
cients  is  given  in  [3 j  and  only  the  applicable  portions  of  the  program 
will  be  mentioned  here. 

The  zonal  mean  height  z  is  defined  as 


S  It 

^  ~  Tirf  1  ° 


(3) 


and  its  zonal  perturbation,  in  accordance  with  equation  (2),  is 

2  -I  (4) 

Expanding  z  in  a  Fourier  series  one  obtains  at  a  given  latitude 

•i 

2  ~  2-  ^  11  [A.  1  l3n  C'c  5  0  (5) 

n  - 1 

where  A  and  B  are  the  coefficients  of  the  nth  harmonic  functions  in  the  • 

Yi  r\  — 

Fourier  expansion  of  z,  and  M  is  the  total  number  of  harmonic  waves  neces¬ 
sary  to  resolve  the  contour  field.  Replacing  and  in  (5)  by  an  am¬ 
plitude  factor  C  ^  and  a  phase  angle  0  ,  (5)  may  be  revrritten  as 
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Equation  (6)  has  been  programmed  for  the  CDC-1604  digital  computer  by 
FNUF,  and  the  program  includes:  a  contour  grid  map  of  the  zonal  pertur¬ 
bation  field  z*,  and  a  tabular  list  of  the  Fourier  parameters  Cn,  ©n  and 
z  for  each  of  27  latitude  circles.  The  last  three  parameters  are  actually 
computed  in  relation  to  a  zonal  band  of  width  equal  to  the  grid  mesh  d  of 
the  FNUF  map  grid,  but  may  be  considered  applicable  at  the  center  of  the 
band . 

Using  spherical  coordinates  h  and  (b  (longitude  and  latitude,  respec¬ 
tively)  the  derivative  operators 
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apply,  and  we  have 
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If  nov/  z  is  replaced  in  (7)  by  its  expansion  in  a  Fourier  series  as  given 
in  equation  (2)  we  obtain 

h  ~  of'  f  c L  ' 111  (  '  "  0t;V 
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Remembering  that  C  ^  and  p  are  independent  of  longitude  4  but  not  of 
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latitude  .f  ,  we  may  perform  the  indicated  differentiations  and  combine 
terms,  with  the  result  that  jf  in  spherical  coordinates  becomes 

h  -  h  i  - j  £l*V«(nA-%)J  ■■%.)(  J"  -  )j 
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Equation  (9)  may  be  shown  to  have  the  form 
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equation  (l( )  may  be  approximated  by 
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Here  the  superscript  bar  denotes  an  average  with  respect  to  n  over  all 
significant  waves  n  =  1 ,  . ,  M. 

The  FNWF  maps  are  stereograp’nic  naps  true  at  latitude  60N,  and  else¬ 
where  the  grid  mesh,  d  =  3&1  km,  is  related  to  the  true  distance  dT  by  the 
relationship  d  =  md_,  where  m  -  (l  +  sin  60)/(l  +  sin,|>)  and  is  the  so- 
called  map  factor. 

The  nartial  derivatives  in  (ll  )  nay  be  evaluated  by  centered  finite 
differences.  The  numerical  values  of  C n  and  were  simply  token  from 
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the  tabular  list  of  computer  output  of  the  harmonic  analysis  at  the  vari¬ 
ous  latitude  circles.  For  example,  d  /  i ' may  be  written  using 
finite  differences  between  adjacent  latitude  circles  (for  which  a  a  cj-,  1  An  ) 
in  the  form 
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The  latitude  circles  of  the  Fourier  analysis  employed  in  the  compita- 
tion  of  J<j  corresponded  to  N  =  12,  15,  and  18  grid  mesh  distances  from  the 
pole,  which  in  turn  corresponded  to  the  latitude  circles:  48. 9H,  38. 6N, 
and  28, 3N  (see  fig,  l).  However,  in  the  computation  of  first  and  second 
derivatives  of  the  nth  Fourier  wave  parameters,  these  parameters  at  lati¬ 
tude  circles  I'M  and  N+1  are  also  needed. 

The  term  vL/  f\  is  more  tractible  in  the  present  case  if  the  finite- 
difference  substitution  ri  A  -  2. d/m  is  made;  thus,  the  centered  deriva¬ 
tive  A  in  finite  differences,  assumes  the  form 


where  T  -  c/  c--:  f  is  the  radius  of  the  latitude  circle,  and  the  subscript 
i+1  is  one  grid  mesh  eastward  of  the  ith  point  alopg  the  latitude  circle. 

Including  the  form  of  ■/■'/.)/.  given  in  (12),  equation  (ll  )  may  now 
be  written  in  the  more  compact  form 


where 
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Fig.  1.  Grid  overlay  depicting  latitude  circles  N=12,  15,  IS  and  the  grid  mesh  length  d.  Geostrophic  relative 
vorticity  is  conjaited  at  the  intersections  of  the  latitude  circles  with  all  5 'longitude  lines  lying 
within  the  area  bounded  by  the  double  -line  meridians. 


The  averaging  with  respect  to  n  in  the  expressions  for  W  and  is 
assumed  to  be  valid,  according  to  the  mean  value  theorem.  Actually  how¬ 
ever  the  average  was  approximated  by  a  simple  arithmetic  mean  of  the  con¬ 
tributions  for  n  =  1 ,  . M.  Thus  it  was  necessary  to  compute  contri¬ 
butions  for  each  n  and  divide  the  result  by  the  number  of  waves. 

Van  Mieghem  |^2  j  points  out  that  waves  of  space  number  higher  than  12 
to  17  contribute  little  to  the  variance  of  the  contour-height  field.  Con¬ 
sequently,  the  Fourier  analysis  was  truncated  at  M  =  12  at  the  two  higher 
latitude  circles,1  although  wave  numbers  1  through  15  were  employed  at  the 
lowest  latitude. 

The  coefficients  and  W2  are  constants  for  a  particular  time  and 
latitude  circle,  and  can  be  computed  directly  from  the  tabulated  harmonic 
data.  These  values  are  listed  in  table  1 .  However,  in  applying  equation 
(13)  it  was  necessary  to  interpolate  values  of  z""  from  the  pert”.rbation 
map  at  the  points  indicated  in  fig,  1  along  the  specified  latitude  circles 
These  points  were  arbitrarily  chosen  at  integral  multiples  of  five  degrees 
of  longitude,  encompassing  the  west-east  extent  of  the  United  States, 
Therefore,  for  ease  in  selecting  data  points,  a  grid  overlay  was  con¬ 
structed  presenting  all  latitude  circles  N?=11 , » . ,  ,,19,  drawn  so  as  to  in¬ 
corporate  the  map  factor. 


At  the  two  higher  latitudes,  H,  was  negative  ensuring  that  lows 
corresponded  to  positive  relative  vorticity.  At  the  lowest  latitude 
(N=18),  it  was  necessary  to  use  M=1 5  in  order  to  ensure  that  this  condi¬ 
tion  was  met. 
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AID  values 


TAB  IE  1  :  Values  of  U,  and  Wx  to  be  applied  in  equation  ( ,  . 

. 1  _  i 

are  in  units  of  10  sec  ft  . 


1“  . .  ""  LATITUDE  CIRCIE 

I - 1 - 


DATE 

1 -  n 

N  =  12 

1 

3 

II 

-1  1 

V/l  j 

1 

1 

r 

N  =  18 

W  =  -.1467 

W 

=  -.1635 

W 

=  -*2696 

3  Nov.  1961 

w  =  -,0079 

1 

•  u 

1 

=  -.0370 

w 

=  -.0303 

W  =  -*1490 

w 

=  -.0617 

w" 

=  -,0420 

4  Nov.  1961 

i 

w  =  -.0311 

w 

=  .0316 

w 

=  .0138 

w  =  -,1355 

w 

=  -*1013 

w 

=  -.1819 

5  Nov.  1961 

11  =  -.0398 

w 

=  -.0230 

•  W 

=  .1682 

1 

i 

j 

| 

j 

»l 

ji 

oj 

1 

1 

i 

w 

1 

"i 

1  ! 

* 

£ 

O! 

°! 

1 

I 

u 

=  -.3510 

6  Nov%  1961 

? 

w  =  .0757 

L  

w 

=  .0349 

w 

=  -.7610 

Mote  that  the  term  (  2  u..,  -  t  V--  in  both  (ll  )  and  (l 3 )  also  varies 

from  point  to  point  around  any  latitude  circle*  However,  the  variation  of 
this  term  was  so  small  that  the  mean  zonal  wind  for  the  latitude  was  used . 
The  mean  zonal  ’wind  was  computed  using  the  latitudinal  gradient  of 
zonal  mean  height  z,  that  is  from 

~  ^  ^  ^ '  *  /  t *  _  1  rfjJ.  -j 

3  r  zd  ’  ; 

again  using  finite  differences* 

Along  each  of  the  three  latitude  circles,  17  values  of  W,  z*  and 
W  (z*( )  were  obtained  for  each  of  the  four  days,  giving  a  total  of 
204  separate  values  of  these  quantities, 

3.  Statistical  results. 

For  precisely  the  same  set  of  ooints  and  days  as  described  above, 
corresponding  values  of  */ere  obtained  from  grid  maps  of  geostrophic 
relative  vorticity  made  available  by  FNTJF .  Their  maps  essentially  pre¬ 
sented  the  field  of  J j  computed  using  their  square-grid  mesh,  according 
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to  the  usual  formula 


I 

f  NWh 


(14) 


where  z  is  the  four-point  space-mean  contour  height  and  d  =  381  km  at 
latitude  6 ON* 

Since  the  averaging  process  in  equation  (ll  )  is  somewhat  approximate, 
it  was  felt  that  better  verification  would  result  if  the  relationship  be¬ 
tween  J  r  and  U  z*  and  W„  ( z ?'  -z  f  , )  was  written 

f  -  rv  *y  r  t  l  a,  a  /  L  i 


J, 


ritHF 


t) 


A.  t- 


)1  (15) 


Here  A^ ,  A. ,  and  A  g  are  regression  coefficients^  to  be  determined  by  a 
least  squares  technique  applied  to  the  204  simultaneous  values  of  W,  zf 
and  W-  (z*(-z5'(  ).  Because  of  occasional  large  gradients  of  '«rNVVf  over  5° 
longitude  intervals,  the  values  of  /  were  smoothed  by  the  three-point 
binomial  running  mean  formula 


*  'h  +  *<>,)  /"/ 


where  is  any  value  of  j(  interpolated  from  the  FNUF  vorticity-grid 
map  at  each  point.  The  regression  equation  actually  involves  the  correla¬ 
tion  of  the  smoothed  J  ^ ..  values  with  the  independent  variables  already 
noted . 

Using  the  BIMD  06  statistical  program  as  adapted  to  the  CDC-1604 
computer,  it  was  necessary  to  type  204  IBM  cards  containing  the  smoothed 
J  ^  values  and  the  simultaneous  values  of  the  independent  variables  of 
the  right  side  of  equation  (15)-  Best-fitting  values  of  the  regression 
coefficients  were 


These  regression  coefficients  are  non-dimensional  since  the  units 
of  U.  zf  and  VI  (z*  -zq. )  are  in  sec-”1. 
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Ac  =  -8.1719, 


A1  =  0.4585, 


and  A?  =  0.0624 


The  standard  error  of  esti.mate"'  and  the  multiple  linear  correlation  co¬ 
efficient  for  the  regression  equation  were  20.3523  and  0.7632  respectively. 

Panofsky  and  Brier  [4,  p.  113]  present  an  F  test  for  the  analysis  of 
the  variance  explained  by  the  regression  equation.  The  F  value  may  be  de¬ 
fined  by  the  ratio  of  the  explained  variance  to  the  residual  variance 
relative  to  the  regression  equation.  This  formula  is 


ft* 

(i.  rx'1)  f> 


(16) 


With  p  =  2,  the  number  of  independent  variables,  and  6}  =  204,  the  number 
of  random  samples,  the  value  o.f  F  from  equation  (l6)  turns  out  to  be 
139.33.  However ,  assuming  normally  distributed  independent  and  dependent 
variables,  the  critical  F  at  the  0.99  significance  level  with  2  and  201 
degrees  of  freedom  is  4.71.  Thus  the  regression  equation  (15)  is  signi¬ 
ficant  at  a  very  high  confidence  level. 

The  partial  correlation  coefficients  of  U  z?  and  W„  (z."  -2*  )  were 
0.7620  and  0.0382,  respectively,  indicating  very  little  contribution  to 
the  regression  equation  for  the  latter  variable.  An  intuitive  explana¬ 
tion  for  the  insignificant  contribution  of  W^  (z*  -z^  )  may  be  obtained 
from  fig.  2,  The  differences,  z?':'.  -z?f  ,  taken  over  the  pressure  centers 
would  be  zero,  whereas  in  actuality  the  largest  absolute  values  of  J ^ 
would  be  found  at  these  points.  likewise  for  a  sinusoidal  contour  wave 
pattern  maximum  magnitudes  of  2‘::'  -2*  would  exist  near  contour  inflection 

4*i  A.-*  * 

points,  where  i  }  -  . 

The  possibility  of  "automatic"  correlation,  viz,  correlation  between 
the  two  independent  variables,  was  examined  and  none  was  found  to  exist. 


All  values  of  are  in  units  x  10“^  sec--' 
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Fig.  2.  Distribution  of  along  a  latitude  circle  between  contour  maxi¬ 
mum  and  minimum,  in  the  case  of  a  sinusoidal  contour  wave. 

The  individual  latitude  circles  were  statistically  examined  in  a  like 
manner.  The  two  northernmost  latitude  circles  gave  a  multiple  linear  cor¬ 
relation  coefficient  of  0.84,  whereas  there  is  no  significant  coi'relation 
between  J  and  the  two  independent  variables  W.zf  and  W,  (z*  -z*  )  at 
latitude  28 ,3N.  This  inconsistency  at  low  latitudes  is  attributable  to 
"noise"  in  the  Fourier  representation  of  the  contour  field.  The  "noise" 
is  due  to  the  small  variance  of  the  height  field  at  this  latitude,  coupled 
with  the  relatively  large  RMS  error  possible  in  the  500  mb  analysis  in  the 
subtropics  during  the  period  of  study.  The  use  of  second  derivatives  of 
the  Fourier  parameters  also  increases  the  amount  of  "noise"  generated  in 
the  determination  of  1  ,  . 

4.  Conclusions. 

Theoretically  there  should  have  been  a  large  positive  multiple  linear 
correlation  between  j  and  W  z V  and  W.  (z;*  -z. ,),  This  proved  to  be  the 

r  i  a.  i  »-/  x 

case  at  latitudes  48. 9N  and  3S.6N  where  the  contour  field  was  well  resolved 
by  twelve  harmonic  waves.  Fowever,  at  the  lowest  latitude  circle,  28. 3F, 
the  several  sources  of  errors  listed  at  the  end  of  section  3  could,  and 
probably  did,  influence  the  low  correlations  that  were  obtained.  These 
were  basically  due  to  the  lack  of  resolution  of  the  contour  field  by  the 
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Fourier  analysis.  However,  the  cane  reasons  would  cause  the  square-grid 
commutation  of.)  >  to  be  subject  to  large  errors. 

Based  upon  the  results  at  the  two  higher  latitudes,  it  seems  apparent 
that  the  use  of  a  Fourier  analysis  of  a  meteorological  field  presents  a 
very  valuable  tool  for  atmospheric  research.  With  the  advent  of  the 
electronic  computer  and  machine -analyzed  contour  charts,  the  use  of  har¬ 
monic  analysis  is  an  effective  and  economically  feasible  approach  to  the 
computation  of  synoptic  aids  for  the  practising  meteorologist. 
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